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Abstract 

In recent years research on credit risk modelling has mainly focused on default probabili- 
ties. Recovery rates are usually modelled independently, quite often they are even assumed 
constant. Then, however, the structural connection between recovery rates and default prob- 
abilities is lost and the tails of the loss distribution can be underestimated considerably. The 
problem of underestimating tail losses becomes even more severe, when calibration issues 
are taken into account. To demonstrate this we choose a Merton-type structural model as 
our reference system. Diffusion and jump-diffusion are considered as underlying processes. 
We run Monte Carlo simulations of this model and calibrate different recovery models to 
the simulation data. For simplicity, we take the default probabilities directly from the sim- 
ulation data. We compare a reduced-form model for recoveries with a constant recovery 
approach. In addition, we consider a functional dependence between recovery rates and 
default probabilities. This dependence can be derived analytically for the diffusion case. We 
find that the constant recovery approach drastically and systematically underestimates the 
tail of the loss distribution. The reduced-form recovery model shows better results, when all 
simulation data is used for calibration. However, if we restrict the simulation data used for 
calibration, the results for the reduced-form model deteriorate. We find the most reliable 
and stable results, when we make use of the functional dependence between recovery rates 
and default probabilities. 
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1. Introduction 



To correctly assess and control exposure to credit risk is a vital question for financial 
institutions. It is also a crucial problem for banking regulation. There are two conceptually 
different approaches to credit risk modelling: structural and r educed- f orm ap proaches. The 
structural models go back to iBlack and Scholesl (119731 ) and iMertonl (119741 ) . The Merton 
model is based on the assumption that a company has a certain amount of zero-coupon 
debt which becomes due at a fixed maturity date. The market value of the company is 
modelled as a stochastic process. A possible default and the associated recovery rate are 
determined directly from this market value at maturity. In the reduced-form approach de- 
fault probabilities and recovery rates are described independently by stochastic models. The 
aim is to describe the dependence of these quantities on com mon covariates or risk f actors . 
For some wel l know n reduced-form model approa c hes see, e.g..lJarrow and Turnbull (1995), 
Jarrow et all (ll997h . iDuffie and Singleton] (Il999h . iHull and White! (l2000h and ISchonbucher 



(120031 ). First Passage Models constitute a third approach which is usually regarded as 
structural, but is better describ ed as a mixed o r pseu do-structural approach. First Passage 
Models were first introduced bv lBlack and Coxi (119761 ). As in the Merton model, the market 
value of a company is modelled as a stochastic process. Default occurs as soon as the market 
value falls below a certain threshold. In contrast to the Merton model, default can occur at 
any time. In this approach, the recovery rate is not determined by the underlying process 
for the market value. Instead, reco yery rates are mode ll ed ind ependently, for example, by 
a reduced-form approach (see e.g., Asvanunt and Staal ( 2009a , bhl. In some cases recovery 



Gieseckel (I2004T )/ 



rates are even assumed constant, for instance, in 

The independent modelling of default and recovery rates can lead to a serious underesti- 
mation of large losses. This situation is particularly troublesome when considering calibra- 
tion to limited default and recovery data. One way to address the scarcity issue of default 
and recovery data is to calibrate to market data of credit contracts and derivatives instead. 
This allows for a self-consistent description of market prices, but may misrepresent the ac- 
tual credit risk. As witnessed in the recent financial crisis, the market as a whole may be 
completely wrong in its pricing consensus. Therefore it is crucial for the risk assessment of 
a credit portfolio to take into account historical data of defaults and recoveries. 

The aim of our study is to discuss general calibration issues that arise from independent 
modelling of default and recovery rates. To this end we use a structural model as our ref- 
erence system. While the original Merton model contains only a diffusion process, we also 
consider a jump-diffusion process for the evolution of the mark et val u e. Fo r related works 
on jum p-d iffusion processes in credi t risk modelling see, e.g., IZhoul ( 120011 ) . ISchafer et al. 



(120071 ) and iKiesel and Schererl (120071 ). We run Monte Carlo simulations of this model and 
calibrate different recovery models to the simulation data. For simplicity, we take the de- 
fault probabilities directly from the simulation data. We compare a reduced-form model 
for recoveries with a constant recovery approach. In addition we consider a functional de- 
pendence between recovery and default rates. This dependence can be derived analytically 
for the diffusion case. We find that the constant recovery approach drastically and system- 
atically underestimates the tail of the loss distribution. The reduced-form recovery model 
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shows better results, when all simulation data is used for calibration. However, if we restrict 
the calibration data, the results for the reduced-form model deteriorate. We find the most 
reliable and stable results, when we make use of the functional dependence between recovery 
and default rates. 

The paper is organized as follows. In Section [2] we discuss the structural model used 
as our reference system. In Section [3] we describe the recovery rate models which will be 
calibrated to the Monte Carlo data. Issues that arise in model calibration are discussed 
in Section HI In Section [5] we compare the different recovery models with respect to three 
central questions: To what extend do the models preserve the dependence of recovery and 
default rates? How robust are the models with respect to calibration issues? And do they 
provide accurate estimates for the risk of large portfolio losses? We summarize our findings 
and conclude in Section El 



2. Structural model as reference system 

In the Merton model, defaults and recoveries are directly determined by an underlying 
market value at maturity. Hence, stochastic modelling of the market value V&(t) of a com- 
pany allows to assess its credit risk. The model assumes that a company k has a certain 
amount of zero-coupon debt with a face value F k . The debt will become due at maturity 
time T. The company defaults if the value of its assets at time T is less than the face value, 
i.e., if Vfc(T) < F k . The recovery rate is then R k = V k (T)/F k and the normalized loss given 
default is 

L\ = 1 - R t = Fk - ' UT) . (1) 
fk 

We denote the loss given default with an asterisk to distinguish it from the loss including 
non-default events. The individual loss can be expressed as 

where G is the Heaviside function. 
A. Model setup 

We model the time evolution of the market value of a single company k by a stochastic 
differential equation of the form 

dV k 



Vk 



fidt + y/c adW m + dJ m + y/l - c adW k + dJ k . (3) 



This is a correlated jump-diffusion process with a deterministic term fidt and a linearly 
correlated diffusion. The Wiener processes dWk and dW m describe the idiosyncratic and the 
market fluctuations, respectively. For simplicity, we choose the drift /x, the volatility a and 
the correlation coefficient c as constant parameters, which are the same for all companies. 
The jump terms are not contained in Merton's original model; they ensure that the default 
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probability does not go to zero as the time to maturity becomes very short. We include two 
jump terms in the stochastic pr ocess, dJk for i diosy ncratic jumps and dJ m for jumps which 



affect the entire market. As in ISchafer et al.l (120071 ) . the jumps are modelled by a Poisson 



process with intensity A. We recall that in such a process the probability function for the 
event to occur n times between zero and the time t is given by 



p Poisson (t)= (ML exp{ _ Xt) _ (4) 



IV. 



The size A of the jump, measured in units of the current asset value V(t), is a random 
variable with a distribution which we have to specify. Jumps can be positive or negative. The 
largest possible negative jump is 100% of the current asset value. Based on this information, 
a possible distribution of the jump size A is a shifted lognormal distribution, A + 1 ~ 
LN(/zj, <tj), with mean \ij and standard deviation oj. 

Without the jump term, the distribution of the asset price Vk(t) is log-normal. The 
lumps render the tails of the asset price distribution fatter. Fat tails are empirically ob- 



served iMailte^n^^Jii^tanieyJ ( 120001 ) . As this clearly affects the loss distribution, we find it 
important to include such jumps. The parameters of the jump process can be adjusted in 
order to match the tail behavior of a given empirical time series of the asset value. We use 
the same parameters for idiosyncratic and market wide jumps. 

B. Monte Carlo simulation 

In the Monte Carlo simulations we consider the stochastic process in Equation fl3]) for 
discrete time increments At = T/N, where the time to maturity T is divided into N steps. 
The market value of company k at maturity is then given by 



iV 

V k (T) = V k {0) Y[ (l + /i At + y/har^VM + dJ m>t + Vl-cae k)t VM + dJ M 
t=i 



(5) 



The random variables r] mjt and Ek,t are independent and are drawn from a normal distribu- 
tion. We consider a homogeneous portfolio of size K with the same parameters for each 
asset process, and with the same face value, Fk = F, and initial market value, Vfc(0) = V$. 
The simulation is run with an inner loop and an outer loop. In the inner loop we simulate 
K = 500 different realizations of £k,t and dJk,t for a single realization of the market terms 
rj m> t and dJ m> t with t — 1, . . . , N. The inner loop can be interpreted as a homogeneous port- 
folio of size K, or simply as an average over the idiosyncratic part of the process. In each 
run of the inner loop, we calculate the market return X m , the number of defaults N-o(X m ) 
and the expected recovery rate (R(X m )). The market return X m is defined as the average 
return at maturity, 

1 K 

X„ 



1 (Vt(T) \ 



For sufficiently large K the idiosyncratic terms average out and the market return X m is 
solely defined by the realizations of 7] myt and dJ m j- This is why we use the market return as 
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a parameter for the other observables. The number of defaults Nv(X m ) simply counts how 
many times the condition Vk(T) < F is fulfilled. We can estimate the default probability as 



P D (X m ) « N D (X m )/K . 



(7) 



We obtain the portfolio loss as the average of individual losses in Equation (j2J), 



(L(X m )> = if^L Jfc . 




(9) 



(R(X m )) = 1 



(jW) 

PD(X m ) 



1 - 



K(L(X m )) 
N B (X m ) 



(10) 



Here, we assume that the number of defaults is strictly non-zero, which is justified for large 
K. The outer loop runs over 10 6 realizations of the market terms, where we obtain different 
values for the market return X m and, consequently, the number of defaults N D (X m ) and the 
expected recovery rate (R(X m )). 

The Monte Carlo simulations are carried out for diffusion and jump diffusion. In this 
paper we do not discuss the parameter dependence of the models or aim at calibrating them 
to a given portfolio. Instead, we only present the results for a single set of parameters 
with economically sensible values. This suffices to make our main argument about different 
recovery model approaches. As correlation coefficient we choose c = 0.5. The parameters 
for the diffusion process are /i = 0.05 and o = 0.15. The additional parameters for the jump 
terms are A = 0.005, /ij = 0.4 and <Tj = 0.3. The initial market value is set to Vq = 100, the 
face value of the zero-coupon bonds is F = 75 with maturity time T = 1. 

3. Recovery rate models 

In our study we want to isolate the influence of the different recovery models on various 
risk measures. Therefore we will not employ a stochastic description of default probabilities. 
Instead we consider the default data from the simulations directly. We compare three 
different recovery rate models: a reduced-form model, constant recoveries and structural 
recoveries. 

A. Reduced-form approach 

Reduced-form models are stochastic models for default probabilities and recovery rates. 
Their aim is to describe the dependence of these quantities on common (macroeconomic) 
covariates or risk factors. In our structural reference model we use a one-factor model with 
a constant correlation of all companies to the market. Therefore we use the market return, 
i.e., the average return of all assets at maturity, as a covariate. 
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Default events are commonly modelled by a doubly stochastic Poisson process or Cox 
process. As mentioned above, we will not consider a stochastic description of default prob- 
abilities in this paper. Instead we consider the default data from the simulations directly. 

Reduced-form recovery rate models describe a monotonous dependence on covariates. 
In our case, for instance, we expect a lower recovery rate for lower market returns. Since 
the natural domain for the recovery rate is the interval [0, 1], a generalized linear model is 
typically used. In this approach a linear dependence on the covariates is transformed by a 
link function to the desired domain. Recovery rate mode ls commonly use eit her a logit or a 



probit link functio n. A logit link function i s used , eg, in ISchonbucherl (120031). a probit lin k 



function is, eg, in lAndersen and Sideniud (120051 ). As pointed out by IChava et al.l (120081 ) 



both link functions lead to very similar results. In our study we focus on the probit model; 
then the reduced-form recovery rate model reads 



(R{X m )) = ${-iX m -6) 



(11) 



with the parameters 7 and S. The qualitative insights presented in this study do not depend 
on the choice of the link function. 



B. Constant recovery rate 

In the simplest reduced-form ansatz the recovery rate is assumed constant, 

(R(X m )) = (R) = const . 



12) 



Although it is obvious that this crude model neglects important structural information and 
therefore cannot be accurate, it is still widely used in research and in practice. The only 
valid statements to make with a constant recovery rate involve setting (Rk) to zero and 
thereby estimating worst case losses. 

C. Structural recovery rate 

In addition to the reduced-form approach, we consider a functional dependence of re- 
covery rates and default probabilities, 



(R(Pd)) = -^exp (-BQ-^Pd) + ^) $ (Q-^Pd) ~ B) . 



(13) 



This result has been derived by lSchafer and Koivusalol (120111 ) in the framework of the Merton 
model for a correlated diffusion process, but it is also a very good approximation for other 
processes. We call this the structural recovery rate model. The relation depends only 
on a single parameter B. In the case of the correlated diffusion, B is determined by the 
parameters of this process. 
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4. Model calibration 



We will use the data obtained in our Monte Carlo simulations for model calibration. 
Since we do not employ a model for default probabilities, these will be directly calculated 
from the number of defaults found in the simulation, 

P D (X m ) « . (14) 

The three different recovery models will be calibrated to the average recovery rates (R(X m )) 
found in the simulation. 

A. Reduced-form recovery rate 

We apply $~ 1 to both sides of Equation ( fTTj) and obtain 

b(X m ) = ((R(X m ))) = - 7 X m - 5 . (15) 

According to our reduced-form model we expect b(X m ) to show a linear dependence. In 
Figure [1] we present the Monte Carlo data of b(X m ) for diffusion and jump-diffusion. In 
addition to the raw data, we also show average values for small intervals of X m . The 
dependence on X m is well defined and rather linear for large negative market returns. This 
is true both for diffusion and for jump-diffusion. For market returns greater than zero, 
the dependence on X m is only visible on average. In single runs of the simulation, the 
recovery rates, and thus also b(X m ), can deviate considerably from the average values. 
The reason for this is rather obvious: positive market returns make defaults less likely 
for all individual companies. Thus, defaults and recoveries are mostly determined by the 
idiosyncratic parts of the processes. The model assumption of a linear behavior of b(X m ) 
is well justified for large negative market returns X m . However, we observe deviations 
of this linear behavior starting already around X m = —0.2. These deviations are most 
visible in the jump-diffusion case, where we even find a negative slope for X m > 0. From 
these findings we can already deduce a qualitative picture of the calibration difficulties that 
may arise. The overall behavior of recovery rates is not captured by the reduced-form 
model. If we calibrate only to data for large negative market returns, we will overestimate 
recovery rates for moderate market returns. And if we only have data available for calibration 
which corresponds to rather moderate market returns, we would overestimate recovery rates 
for large negative market returns. In both cases, we overestimate recovery rates in some 
situations, and thus underestimate the corresponding losses. Here we choose as an upper 
threshold for X m and study the calibration results in dependence of a lower threshold. The 
motivation for this is the typical abundance of data for moderate market returns, while 
large negative market returns occur only scarcely. We calibrate to the local average values 
in order to give equal weights to different values of X m in the regression. 

Another problem becomes obvious when we compare the diffusion and the jump-diffusion 
case: The quality of the reduced-form recovery model depends on the underlying process. 
In other words, a model that works well for one set of credit contracts might not work so 
well for another. Non-stationarity in financial time series may also cause the quality of a 
reduced-form recovery model to degrade over time. 
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Figure 1: Regression data b(X m ) — ((R(X m ))) for the reduced-form recovery model. The model 
assumes a linear dependence, which is justified for large negative market returns X m . Results are shown for 
diffusion (left), and for jump-diffusion (right). 



B. Constant recovery rate 

The constant recovery rate is determined as the average value of (R(X m )) for a given 
range of market returns X m . The reasons for choosing an upper and a lower threshold on 
X m have been discussed above. The upper threshold X m < addresses the bias issue of the 
estimation at least to some degree. 

C. Structural recovery rate 

In our structural recovery approach, we do not fit a covariate dependence of the recovery 
rate explicitly. Instead, we focus on how recovery rates depend on default probabilities. 
The functional dependence of recovery rate on default probability depends only on a single 
parameter B. 

(R k (P D )) = ^-ex V (^-B^-\P D ) + ^B 2 ^(^- 1 (P D )-B) . (16) 

In the Merton model for a correlated diffusion process the parameter B is given by the 
parameters of this process, 

B=y/(1- c)a*T . (17) 

We use a least square method to estimate B from the simulation data for Po(X m ) and 
(Rk{X m )). Again, the data are filtered with respect to an upper and a lower threshold on the 
m arket return X m . The converge nce of the regression for B is strong. As already pointed out 



in 



Schafer and Koivusalol (120111 ). this functional dependence does not critically depend on 
the underlying process. In fact, it has been shown to work well for diffusion, jump diffusion 
and GARCH(1,1). Therefore we do not expect a critical dependence of the calibration results 
on the range for X m . We can further improve the stability of the calibration if we regress 



S 



the data of average losses directly. Inserting the functional dependence (TT3|) into Equation 
CO!) yields 

(L(P D ))=P D -exp^- J B$- 1 (P D ) + i J B 2 ^$($- 1 (P D )- J B) . (18) 

Now the strongly fluctuating results for small default probabilities are suppressed in the 
regression, see Figure [2J This leads to an improvement in calibration stability. 
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Figure 2: Dependence of average losses (L(X m )} on default probabilities Pu(X m ). Results are shown for 
diffusion (left), and for jump-diffusion (right). 



5. Comparison of model results 

A. Dependence of recovery rates and default probabilities 

Let us first turn our attention to the dependence between recovery rates and default 
probabilities. Figure [3] shows scatter plots of the Monte Carlo results for PD(X m ) and 
(R(X m )). We compare these with the results for the three different recovery rate models, 
where P^(X m ) is again taken directly from the simulation. Here, the model parameters for 
the recovery rate are calibrated to the data with X m < 0. It is obvious that the assumption 
of a constant recovery rate is not justified at all. Moreover, the estimation of this constant is 
biased due to the relative abundance of data for small default probability and high recovery 
rate. Although this bias could be taken into account, the results for the constant recovery 
rate model appear arbitrary and neglect important structural information. For the reduced- 
form recovery, we notice that the interdependence of recovery and default rates is preserved 
at least qualitatively. We observe larger deviations for the jump-diffusion. This is to be 
expected, since the recovery model does not describe the simulation results over the whole 
data range, see our discussion in the previous section. The qualitative agreement of the 
reduced-form model results and the simulation data can be attributed to the fact that 
both P D (X m ) and (R(X m )) describe the dependence on the same covariate. This is in 
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accordance with findings of IChava et al.l ( 120081 ) who pointed out that reduced-form models 
can yield considerably better results, if defaults and recoveries are modelled and calibrated 
with respect to the same set of covariates. 

The structural recovery describes the data very well over the whole data range. This is 
not surprising for the diffusion process, after all the functional dependence has been derived 
for this case. It is noteworthy, however, that we find such a perfect agreement for t he jum p- 
diffusion as well. This result was previously mentioned in ISchafer and Koivusalol (j201lh . 
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Figure 3: Dependence of recovery rates (R(X m )) on default probability PD(^ m )- The Monte Carlo results 
are compared to structural (solid line), reduced-form (dashed line) and constant recovery rate (dashed- 
dotted line). Results are shown for diffusion (left), and for jump-diffusion (right). 



B. Value at Risk and Expected Tail Loss 

Finally, we want to test the quality of the model results on two risk measures that 
capture the tail behavior of the loss distribution. For this purpose, we consider the Value 
at Risk and the Expected Tail Loss, or Ex pected Shortfa l l. Bot h measures are commonl y 
used in credit risk manan gement, see, e.g., lArtzner et al.l (119971); iFrev a nd McNeil (2002); 



Yamai and Yoshibal (120051 ) . and in financial regulations, see lBasel Committee on Banking Supervision 
2006h . For details on the calculation of these risk measures we refer to |Appendix A| As 



mentioned above, we are considering a situation where the default probabilities are well 
described by a model — or in our case by the simulation itself. Hence, we focus entirely on 
the influence of the different recovery rate models. In particular, we are interested in the 
impact of limited calibration data on the model results. To this end we filter the simulation 
data with respect to the market return X m . As before we keep as an upper threshold, 
because default and recovery rates do not show a clear dependence on the market return for 
positive X m . Additionally, we set a lower threshold on X m and study the results for VaR 
and ETL in dependence of this lower threshold. This corresponds to the relative abundance 
of empirical data for moderate market returns, while large negative market returns occur 
only scarcely. Frankly, we ask to which extend the different models can describe extreme 
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events when they are only calibrated with default and recovery data from a relatively calm 
period. 

Figure H] shows the model results for the Value at Risk VaR a at the confidence level 
a = 0.99. The results are normalized to the corresponding empirical value which is computed 
for the entire data set of the simulation. These empirical values are VaRo.99 = 0.013 for the 
diffusion and VaRo.99 = 0.015 for the jump-diffusion. Not surprisingly, we find the worst 
results for the constant recovery. The Value at Risk is underestimated by more than 30% in 
the diffusion case and by almost 20% in the jump-diffusion case. The better result for the 
jump-diffusion is misleading, however, as we will see later on. The reduced-form recovery 
yields quite reasonable results in the diffusion case. It overestimates the Value at Risk by 
roughly 8% when all data with X m < is taken into account. We see a slight dependence 
on the lower threshold for X m . This starts only at X m f» —0.35, because this was the lowest 
market return observed in the simulation. In the jump-diffusion case, the reduced-form 
model yields very unstable results, ranging from a considerable overestimation to a severe 
underestimation of the Value at Risk. For both processes, the structural recovery model 
provides rather stable results with very little deviation from the empirical values. 

In Figure [5] we present the model results for the Expected Tail Loss ETL a at the confi- 
dence level a = 0.99. Again, we normalized the results to the empirical value for the entire 
data set of the simulation. These empirical values are ETL0.99 = 0.0238 for the diffusion and 
ETL0.99 = 0.0379 for the jump-diffusion. The Expected Tail Los s is better suited to d escribe 
the tail behavior of the loss distibution, as first pointed out by lArtzner et al.l ( 119971 ). Here 
we see that the constant recovery rate does in fact not provide a better description of the 
jump-diffusion than of the diffusion. The results for the Value at Risk are more coincidental, 
they depend critically on the chosen confidence level a. The reduced-form recovery again 
shows rather good results for the diffusion process. However, the Expected Tail Loss reveals 
even more clearly that the reduced-form recovery severly underestimates large losses, if the 
calibration data is restricted. Moreover, we observe a very strong dependence of the results 
on the lower threshold for the market returns. These unstable results are contrasted by the 
stable results for the structural recovery rate. The Expected Tail Loss confirms that this 
model delivers a very good and robust estimate for large portfolio losses. 



6. Conclusions 

The relation between recovery rates and default probabilities has a pronounced effect 
on the tail of the portfolio loss distribution. Yet both quantities are often modelled in- 
dependently in current credit risk models. In this paper we consider a Merton model as 
a reference system with diffusion and jump-diffusion as underlying processes. Three dif- 
ferent recovery models are calibrated to Monte Carlo simulations of this reference system: 
a reduced-form recovery model, the special case of constant recoveries, and the structural 
recovery model. A functional dependence between default and recovery rates can be derived 
for the Merton model with correlated diffusion. The structural recovery model employs this 
same functional dependence regardless of the underlying process. We study to what extend 
the different recovery models preserve the relation between default and recovery rates. Most 
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Figure 4: Value at Risk VaR Q at the confidence level a = 0.99. The results are normalized to the 
Value at Risk computed for the entire data set of the simulations. Thus, values above 1 correspond to an 
overestimation, values below 1 to an underestimation by the model. Results are shown for diffusion (left), 
and for jump-diffusion (right). 
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Figure 5: Expected Tail Loss ETL Q at the confidence level a = 0.99. The results are normalized to the ETL 
computed for the entire data set of the simulations. Thus, values above 1 correspond to an overestimation, 
values below f to an underestimation by the model. Results are shown for diffusion (left), and for jump- 
diffusion (right). 



importantly, we compare how well they reproduce two measures of tail risk, the Value at 
Risk and the Expected Tail Loss. Furthermore, we discuss calibration issues that may arise 
due to a lack of information. To this end we confine the calibration to a certain range of 
the simulation data. A lower threshold on the market return implies that extremely adverse 
market situations are not contained in the historical data sample on defaults and recoveries. 

It comes as no surprise that the constant recovery rate model shows the worst perfor- 
mance, since it does not include a dependence on covariates or default probabilities. In 
addition, the estimate of the average recovery rate from historical or simulation data is bi- 



12 



ased due to the relative abundance of small losses. As a consequence this model shows very- 
poor results for Value at Risk and Expected Tail Loss. Nonetheless we included this model 
in our discussion, as it is still prevalent in practical applications. 

The results for the reduced-form recovery model are considerably better. The generalized 
linear dependence on market returns works reasonably well, especially when considering 
large negative market returns. However, it does not provide a good description over entire 
data range. Furthermore, the quality of the model critically depends on the underlying 
process. The model is able to reproduce the dependence of default and recovery rates at 
least qualitatively. Larger deviations are observed for the jump diffusion process. We find a 
critical dependence of the model results on the data range used for calibration. This leads 
to very unstable estimates for Value at Risk and Expected Tail Loss. 

The structural recovery model accurately describes the dependence of default and recov- 
ery rates for different processes. It has only a single parameter and is calibrated with respect 
to the default rate dependence. The model yields nearly perfect estimates of Value at Risk 
and Expected Tail Loss. These results show almost no dependence on the data range used 
for calibration. 



In accordance with IChava et al.l (120081 ) we conclude that a reduced-form ansatz can 
work well, if the same covariates are used for both default and recover rates. However, an 
additional model risk is involved in the reduced-form recovery with respect to the underlying 
process. The structural recovery avoids this model risk because it does not rely on covariates, 
instead it directly describes the dependence of recovery rates on default probabilities. Thus, 
it can be easily used with any default model. 
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Appendix A. Calculation of risk measures 

In the reduced-form approach default probabilities and recovery rates are modelled in 
dependence of covariates, in our case, in dependence of the market return X m . The average 
loss or loss of a homogeneous portfolio is then 

(L(X m )) = P D (X m ) (1 - (R(X m )}) . (A.l) 

There is an abundance of empirical data on covariates, we can therefore assume that their 
distribution is well understood and modelled. In the case of a correlated diffusion, we know 
that the market return is log- normal distributed. The jump-diffusion process renders the tail 
of the distribution fatter, corresponding more closely to the situation on the stock market. 
When the PDF px m (X m ) is known, we can use the substitution 

\ PL (L)dL\ = \ Px jX m )dX m \ (A.2) 
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to calculate the Value at Risk and the Expected Tail Loss. Here we simplified the notation 
to L instead of (L(X m )). With the cumulative distribution function 



FxM = J P xAX m )dX m (A.3) 

— oo 

we can express the Value at Risk VaR Q as 

VaR a = L(F^(l-a)), (A.4) 

where the function L is the V m -dependence given in Equation (1A.1I) . The value F% (1 — a) 
is the 1 — a quantile of the market return. For numerical data it can be calculated directly. 
The Expected Tail Loss is calculated as 

1 F^Jl-a) 

ETL a = — / Lp L (L)dL = - [ L(X m )p x jX m )dX m . (A.5) 
a J a J 

VaR a — oo 

For numerical data we can determine the Expected Tail Loss as 

ETL a = mean {L{X m ) \ X m < F x x Jl - a)} . (A.6) 

The same line of reasoning applies for the structural recovery approach. The average loss is 
then given by 

(L(P D )) = P D (1 - (R(P D ))) , (A.7) 

with the expected recovery rate according to Equation ffT3]) . When we know the distribution 
of default probabilities, e.g., by means of modelling, by a scenario analysis or from Monte 
Carlo simulations, we can express the Value at Risk as 

VaR Q = L(F^(a)) , (A.8) 

where Fp D is the cumulative distribution function of the default probability. For the Ex- 
pected Tail Loss we find 

i 

ETL Q = — / L(P D )p PD (P D )dP D . (A.9) 
a J 

For numerical data we can evaluate it as 

ETL Q = mean (L(P D ) | P D > F p ^(a)} . (A.10) 
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